Surface Roughness Dominated Pinning Mechanism of Magnetic 
Vortices in Soft Ferromagnetic Films 

T. Y. Chen, M. J. Erickson, and P. A. Crowell 

School of Physics and Astronomy, University of Minnesota, 
116 Church St. SE, Minneapolis, MN 55455 

C. Leighton 

Department of Chemical Engineering and Materials Science, 
University of Minnesota, ^21 Washington Ave. SE, Minneapolis, MN 55455 

Abstract 

Although pinning of domain walls in ferromagnets is ubiquitous, the absence of an appropriate 
characterization tool has limited the ability to correlate the physical and magnetic microstructures 
of ferromagnetic films with specific pinning mechanisms. Here, we show that the pinning of a 
magnetic vortex, the simplest possible domain structure in soft ferromagnets, is strongly correlated 
with surface roughness, and we make a quantitative comparison of the pinning energy and spatial 
range in films of various thickness. The results demonstrate that thickness fluctuations on the 
lateral length scale of the vortex core diameter, i.e. an effective roughness at a specific length 
scale, provides the dominant pinning mechanism. We argue that this mechanism will be important 
in virtually any soft ferromagnetic film. 
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The pinning of domain walls in ferromagnets is attributed to the interactions between 
the domain structure and local fluctuations of magnetic properties due to defects. Possible 
sources of defects in polycrystalline ferromagnets include point defects (e.g. impurities, va- 
cancies, and nonmagnetic inclusions), line defects (e.g. dislocations), surface imperfections 
(e.g. roughness), and random anisotropies. Because multiple types of defects coexist in a 
given material, it is problematic to identify exactly which ones dominate the pinning pro- 
cess. It has been practically impossible to identify individual pinning sites, and therefore 
studies of domain wall pinning have focused on collective effects jl-4|. This approach is 
not adequate for applications of domain wall-based devices in which pinning must be pre- 
cisely engineered. An alternative approach is to study pinning in simple, albeit non-uniform 
magnetic structures, where individual pinning sites can be readily identified. Recent stud- 
ies of single magnetic vortices in ferromagnetic disks provide excellent examples, showing 
discontinuous vortex motion as a function of the applied magnetic field 0, defect-induced 



enhancement of the gyrotropic frequency 



harmonic pinning potentials {9]. It is thus clear from these studies that pinning of a single 
vortex can in fact be probed via vortex dynamics, although the precise pinning mechanism 
has remained largely mysterious. 

In this Letter, we report on the dominant pinning mechanism for a single magnetic 
vortex in soft ferromagnetic permalloy (NigoFe2o) films, quantifying pinning energies and 
spatial ranges as a function of film thickness. We show that the measured pinning range, 
approximately 20 nm, is nearly identical to the vortex core diameter, demonstrating that 
the pinning defects interact only with the core of the vortex. We further show, using the 
thickness dependence of the pinning energy, that the dominant pinning defects are located on 
the surfaces. We demonstrate quantitatively that the pinning is correlated not with the root 
mean square (RMS) of the surface roughness but rather with the roughness on the lateral 
length scale of the core diameter, i.e. an effective roughness. We argue that this vortex- 
pinning mechanism will be important in virtually any soft ferromagnetic film. Our findings 
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Stable vortices are obtained in micron size ferromagnetic disks in which the magnetization 
curls in the plane of the disk, with the exception of the disk center, where the magnetization 
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FIG. 1. (Color online) Contour maps of the gyrotropic frequency fa as a function of the in-plane 
static field for disk thicknesses of (a) 20 nm, (b) 35 nm, (c) 50 nm, (d) 65 nm, (e) 80 nm, (f) 100 
nm, and (g) 130 nm. High frequency areas in the contour maps correspond to pinning sites, as 
indicated by the dashed circle in (c) for example. AH is the radius of a pinning site as defined in 
the text. 



orients out of the plane within a core region on the order of 10 nm in diameter. In our case 
the magnetic disks were patterned from polycrystalline permalloy (NigoFe2o) films, which 
were grown on Si substrates with a SiN buffer layer by dc magnetron sputtering at 100 W 
(0.1 nm/s) in 2.5 mTorr Ar, at ambient temperature. The film thicknesses studied were 
20, 35, 50, 65, 80, 100, and 130 nm. On each permalloy film, a 70 nm thick Ti layer was 
deposited as a hard mask, which was patterned into l-/im-diameter disks using electron 
beam lithography. The disk patterns were then transferred to the NigoFe2o layer by Ar ion 
milling. 

We used time-resolved Kerr microscopy (TRKM) to measure the vortex gyrotropic mode 
of individual magnetic disks, i.e. the lowest frequency excitation of the vortex [l^]. We first 
mapped the spatial distribution of pinning defects by measuring the gyrotropic frequency 
fa as a function of the orthogonal in-plane d.c. magnetic fields js, 7]. The in-plane fields 
were varied in increments of 5 Oe over a range of 100 Oe x 200 Oe, which displaces the 
vortex core over a 110 x 220 nm 2 spatial region around the center of the disk. For each set 
of static magnetic field values, the gyrotropic mode was excited by a magnetic field pulse 
with a temporal width less than 120 ps and an amplitude of 5 Oe oriented in the plane of 
the disk, and the resulting gyrotropic frequency f G was measured. Contour maps of fa as 
a function of the static fields are shown in Fig. [TJ f G is represented by a color scale, and 



pinning sites appear as localized regions of high fa- 

We characterized each pinning site in Fig. [T] via two quantities, the pinned frequency 
f P i n , which is the highest frequency within each point-like area, and the depinning field AH, 
where 2AH is the full width at half maximum (FWHM) of the fg peak, averaged from 
the two orthogonal field directions. The averaged pinning-site characteristics ((f P i n ) and 
(AH)) for each sample are shown in Figs. [2^a) and (b) respectively. It is clear that (fp% n ) is 
significantly higher in thinner disks, while (AH) is relatively insensitive to thickness L. 

Also shown in Fig. |2Ja), is the "unpinned" gyrotropic frequency L, which is obtained by 
measuring fa at higher amplitude to remove the influence of pinning J9| . L increases with L 
as expected from analytical models and micromagnetic simulations |13l-ll6|. It is surprising 
however, that the enhancement of fa due to pinning ((f P i n ) — f u ) varies approximately as 
1/L. As will be discussed below, the enhancement of fo is associated with the lateral range 
of pinning, while the 1/L trend reflects the existence of a single length scale characterizing 
the pinning interaction. 

Within a simple model, the measured AH and f P i n allow us to determine the physical 
properties of each pinning site, including the pinning energy E p i n and the pinning range D p i n . 
This can be done within the "two- vortices" model [14 1 developed to describe vortex dynamics 
in a disk, considering the local pinning potential created by a defect. The core is assumed 
rigid, and the geometric confinement W is approximated by a parabolic potential, W(r) = 
k u r 2 /2, where r is the distance between the core and the disk center, k u = M^L^n / 'xo 
is the unpinned stiffness, M s is the saturation magnetization, xo = R/10L is the vortex 
susceptibility, and £ = 2/3 is a model-dependent constant. The applied in-plane magnetic 
field H changes the potential energy by Hfi(r), where /i(r) = £M s 7irRL, and R is the 
radius of the disk. Thus, in an applied field the core moves to a new equilibrium position 
r(H) = xoRH/(M s £) for an unpinned vortex. In contrast, for a pinned vortex the core is 
trapped by a local pinning potential W p i n (r), unless a sufficiently large magnetic field ~ 
AH is applied to overcome the energy barrier. This is illustrated in the inset of Fig. [2] For 
a pinning site in the middle of the disk, we can estimate the pinning energy E p i n , i.e. the 
energy barrier, from AH using 

= \k u r{AH) 2 . (1) 

This can be generalized to a pinning site at any location in the disk by shifting the origin. 
As shown in Fig. |2] (c), (E pin ) is constant at small thickness and then drifts slightly upward 
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FIG. 2. (Color online) Thickness dependence of (a) the pinned (f p i n ) and the unpinned (f u ) 
gyrotropic frequency, (b) the averaged depinning field AH, (c) the pinning energy E p i n , and (d) 
the pinning range D p i n . The solid line in (d) is the size of core obtained from the micromagnetic 
simulation. The inset shows a schematic of how the total potential energy changes in an applied 
field. 



with L. If the pinning effect were dominated by point defects in the bulk of the film, we 
would expect the pinning energy to scale with the number of defects intercepted by the core 
as it traverses the film. The pinning energy would then be proportional to L, which is not 
observed. Rather, the thickness dependence suggests that the observed vortex pinning is 
dominated by surface defects, motivating a detailed analysis of surface roughness as provided 
below. 

To determine D p i n , we consider the simplest form of a pinning potential, W P i n (r) = 
Akr 2 /2, where Ak = k pin — k u = 2 r n(f pin — f u )G, and G = 2nLM s /'y is the gyroconstant 
We define D pin as the diameter of the pinning potential, so that W P i n (D pin /2) = E pin , 

and 
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As shown in Fig. EJ^d), D p i n approximately matches the diameter of the vortex core, which 
is indicated by the solid curve. The core diameter is obtained from a micromagnetic simu- 
lation using typical material constants for permalloy, including M s = 800 emu/cm 3 and the 
exchange constant A = 1.05 x 10~ 6 ergs/cm. Because the core diameter at the disk surface 
is smaller then in the equatorial plane of the disk [17J, in the simulation we divide the disk 
into 5 layers, and we define the core diameter by the radius of the maximum gradient of 
the z-component magnetization at the surface of the disk. It is likely that local defects pin 
the vortex core region more effectively than the other regions of the vortex due to the large 
energy density within the core, where the magnetization gradient is large. Imperfections on 
length scales similar to the core diameter therefore cause the strongest pinning effects. 
It should be noted that the above explanation for the observed pinning range is consistent 



with the consensus based on models of domain wall pinning 181 ] . Specifically, pinning of 
domain wall motion is known to be most effective at defects with dimensions comparable to 
the wall width. However, to our knowledge this limiting-defect-size effect has been shown 
in experiments only through collective effects in which the highest coercivities 

were observed when the grain or inclusion size in the film matched the estimated domain 
wall width. Here, because we directly identify the spatial range of the interaction between 
a single vortex and an individual pinning site, which is set by the core diameter (Fig. 121(d)), 
our findings represent strong evidence for the predicted limiting-defect-size effect. 

We now turn to analysis of the surface roughness, which was characterized using tapping- 
mode atomic force microscopy (AFM). Representative AFM images of the 20-nm-thick and 
130-nm-thick samples are shown in the inset of Fig. ]3](a). The grain size and roughness are 
found to increase with thickness, as is typical. This can be seen in Fig. E(a), which shows the 
roughness power spectral density (PSD). By integrating the spectra, we can obtain the RMS 
value of the roughness (Jrms- To determine the contribution to the roughness from length 
scales on the order of the core diameter, we integrate the roughness PSD with a weighting 
function, shown as the Gaussian solid curve in Fig. |3]^a). This function is determined 
experimentally from the probability distribution of Tr/D pin , which is shown in Fig. El^b). 
We use 7r/D P i n as the relevant length scale because the vortex core will be accommodated 
most easily when fluctuations in the pinning potential are twice its diameter. Using this 
procedure, we obtain an effective roughness a e ff for each sample, as shown in Fig. M^c). a e ff 
increases from approximately 2 to 4 nm as L increases from 20 to 130 nm. 
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FIG. 3. (Color online) (a) Roughness power spectral density (PSD) as a function of wave-vector 
(lines), obtained from the Fourier transforms of atomic force micrographs (AFM). (b) Probability 
distribution of the pinning range as a function of 7r/D p i n . Data are obtained from the pinning sites 
in Fig.l. This distribution is used to determine the weighting function shown by the solid Gaussian 
curve in (a), (c) Effective roughness a e ff versus the disk thickness L. The two points shown in 
open symbols are additional 50 nm thick films prepared with different roughness characteristics. 



The relationship between pinning and roughness is shown in Figs. H](a) and (b), which 
show E pin as a function of orms and cr e ff, respectively. Besides the samples discussed thus 
far (filled circles), two additional 50 nm thick samples are shown for comparison, indicated 
by open circles for 50 nm sample #2 and open triangles for 50 nm sample #3. Sample #2 
has smaller roughness compared to the original 50 nm sample (#1), and was prepared in a 
different deposition run. Sample #3 has larger grain size (~100 nm) than samples $T and 
#2 (~30 nm), as it was deposited at an elevated substrate temperature of 250 C [7j. As can 
be seen in Fig. S](a), there is no clear linear correlation between (Jrms an d E pin . However, as 
shown in Fig. IH(b), it is clear that E pin scales linearly with <J e ff, and the linear fit, shown as 
a solid line, intercepts with the y-axis approximately at zero. This result indicates that the 
effective roughness <7 e // is the dominant vortex pinning mechanism for all of these permalloy 
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FIG. 4. (Color online) Pinning energy versus (a) the RMS roughness, and (b) the effective roughness 
(on the length scale of the core diameter). The solid line is a linear fit. (c) Pinning energy 
normalized by the effective roughness versus disk thickness. Data obtained on the 50 nm thick 
sample #2 are indicated by open circles. Data obtained on the 50 nm thick sample #3 are indicated 
by open triangles. The dashed line indicates the average value of E p i n /a e ff. 

films. The correlation between E P i n and a e ff also explains the thickness dependence of E P i n , 
shown in Fig. |2]^c). This is due to the fact that a e ff is larger for thicker disks (Fig. L3Ic)), as 
evidenced by the constant value of E pin /a e ff versus thickness, as shown in Fig. LUc). 

We can estimate how small the effective roughness would need to be in order to avoid 
the roughness- induced pinning mechanism. We consider that a pinning site is unimportant 
only if depinning of a vortex from that site can be thermally activated at room temperature 
on a typical laboratory time scale, i.e. 1 sec. With an attempt frequency of 0.5 GHz, 
which is set by the gyrotropic frequency, the critical pinning energy turns out to be 0.5 
eV. Therefore, from the data in Fig. H^b), we determine that the effective roughness would 
need to be smaller than 1.5 nm to avoid vortex pinning. To put this effective roughness 
value in perspective, we consider its implication for films with various growth modes. In 
the Volmer- Weber or island growth mode (often relevant to sputtered polycrystalline metal 
films), roughness is significant and the lateral correlation length is directly linked to the 
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grain size. Given that this length scale in most practical situations is of the same order 
of magnitude as the physically relevant pinning length scale (i.e. 20 nm, the vortex core 
diameter), the constraint cr e ff < 1.5 nm is a very stringent one, requiring grain sizes very 
different from the core diameter. In Frank-van der Merwe (i.e. layer-by-layer) or step- 
flow growth modes (potentially relevant to MBE-grown epitaxial metal films) the lateral 
correlation length of the roughness is set by the mean terrace width, and thus the vicinality 
of the substrate surface. Satisfying cr e // < 1.5 nm may be possible, but even in this case 
it would require specific tailoring of the vicinality and terrace width to avoid the scale of 
the core diameter. Thus, we expect that the surface roughness pinning mechanism plays 
an important role for vortex pinning in virtually all soft ferromagnetic films. Similar length 
scale arguments apply to edge roughness in patterned thin film devices, as the length scale 
associated with the patterning technique, such as electron beam lithography, is again likely 
on a similar scale to the vortex core. 

In conclusion, we have shown that the dynamics of a single vortex allow us to quantify 
both the energies and length scales associated with individual pinning sites. The dominant 
pinning mechanism is the interaction between the vortex core and surface roughness on a 
lateral length scale set by the core diameter. We suggest that this mechanism determines 
the minimum pinning energy for vortex motion in most soft ferromagnetic thin-film devices, 
such as vortex domain walls in nanowires. 
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